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Anomalous critical exponents in the anisotropic Ashkin-Teller model 
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We perform a rigorous computation of the specific heat of the Ashkin-Teller model in the case of 
small interaction and we explain how the universality-nonuniversality crossover is realized when the 
isotropic limit is reached. We prove that, even in the region where universality for the specific heat 
holds, anomalous critical exponents appear: for instance we predict the existence of a previously 
unknown anomalous exponent, continuously varying with the strength of the interaction, describing 
how the difference between the critical temperatures rescales with the anisotropy parameter. 
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More than half a century ago Ashkin and Teller 
introduced their model as a generalization of the Ising 
model to a four component system. It describes a bidi- 
mensional lattice, each site of which is occupied by one 
of four kinds of atoms: A, B, C, D. Two neighbouring 
atoms interact with an energy: eo for AA, BB, CC, DD; 
ei for AB, CD; 62 for AC, BD; and £3 for AD, BC. Fan 
[2] has shown that the AT model can be also written in 



terms of Ising variables — ±.l,a\t' — ±1 located at 
each site of the lattice; its Hamiltonian can be written, 
if X, y are nearest neighbor sites as: 

(1) 

with J(i) = /3(eo + £1 - £2 - £3)74, J(^' = /3(£o + £2 - 
£1 — £3)74, I — /3(£o + £3 — £1 — £2)74 and (3 is the in- 
verse temperature. The AT model is then equivalent to 
two Ising models coupled by an interaction quartic in 
the spins; the case in which the two Ising subsystems are 
identical J^^-' = J^^-' is called isotropic, the opposite case 
anisotropic. When the coupling Z is = 0, Ashkin-Teller 
(AT) reduces to two independent Ising models and it has 
of course two critical temperatures if j'^^^ ^ 

Layers of atoms and molecules adsorbed on clean sur- 
faces, like submonolayers of Se adsorbed on Ni, are be- 
lieved to constitute physical realizations of the AT model 
1^ IS i theoretical results on it can explain the phase 
diagrams of such systems, experimentally obtained by 
means of electron diffraction techniques. As for the Ising 
model, the importance of AT is also in providing a con- 
ceptual laboratory in which the higly non trivial phe- 
nomenon of phase transitions can be understood quan- 
titatively in a relatively manegeable model; in particu- 
lar it has attracted great theoretical interest because is 
a simple and non trivial generalization of the Ising and 
four-state Potts models. 

Contrary to many 2d models in statistical mechanics 
like the Ising, the 6 or the 8 vertex models j^, in which 
remarkable exact solutions give us very detailed informa- 
tions about the behaviour of thermodinamical functions. 



.(2) 



there are no exact results on the AT model except for 
the trivial I = case. It is believed that the AT 
has two critical temperatures for J^^^ ^ J^^) which co- 
incide at the isotropic point J^^^ — J^'^\ Moreover it 
was conjectured by Kadanoff and Baxter that the 
critical properties in the anisotropic and in the isotropic 
case are completely different; in the first case the criti- 
cal behaviour should be described in terms of universal 
critical indices (identical to those of the 2d Ising model) 
while in the isotropic case the critical behaviour should 
be nonunivcrsal and described in terms of indexes which 
are non trivial functions of L In other words, the AT 
model should exhibit a universal-nonuniversal crossover 
when the isotropic point is reached. 

Evidence for the validity of nonuniversal behaviour in 
the isotropic case was given in (using second order 
renormalization group arguments) and in |lOl [Tll | (by a 
heuristic mapping into the massive Luttinger model de- 
scribing one dimensional interacting fermions in the con- 
tinuum). The anisotropic case was studied numerically 
by Migdal-Kadanoff Renormalization Group 's'l, Monte 
Carlo Renormalization group 12], finite size scaling [isj : 
such results give evidence of the fact that, far away from 
the isotropic point, AT has two critical points and be- 
longs to the same universality class of the Ising model 
but give essentially no informations on the critical be- 
haviour when the anisotropy is small. 

In this Letter we present a rigorous derivation of the 
specific heat for the AT model, valid for small interac- 
tion I and any anisotropy. We find indeed that in the 
anisotropic case the specific heat is singular in corre- 
spondence of two critical temperatures, and the diver- 
gence is logarithmic as in the Ising model, in agreement 
with universality hypothesis. Nevertheless even in the 
region where universality holds, anomalous critical expo- 
nents appear; for instance the difference between the two 
critical temperatures rescale with the anisotropy param- 
eter with a nonuniversal critical exponent. The presence 
of such critical exponents also in the universality region 
clarify how the universal-nonuniversal crossover is real- 
ized when the isotropic limit is reached. 
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Such results are found by the new methods introduced 
in 01 and 0| to study 2d statistical mechanics mod- 
els which can be considered perturbations of the Ising 
model. These methods take advantage from the fact that 
such systems can be exactly mapped in models of weakly 
interacting relativistic fermions in rf = 1 + 1 on a lat- 
tice. The mapping was known since long time (see Ref. 
0^ Jifij J.^, 18] ); however in recent years a great progress 
has been achieved in the evaluation of Grassmann in- 
tegrals of interacting models, in the context of quantum 
field theory and solid state physics (see Ref. flO, 20, 21j), 
and one can take advantage of this new technology to get 
informations about 2d statistical mechanics models. This 
provides the only method to get rigorous quantitative in- 
formations on the critical properties of such systems if 
an exact solution is lacking, as in the present case. The 
algorithm is based on multiscale analysis and allows us to 
prove convergence of the expansion for the energy-energy 
correlation functions and for the specific heat up to the 
critical temperature; essential ingredients of our analysis 
are cancellations due to anticommutativity of fermionic 
variables and approximate Ward identities p^ . guaran- 
teeing that the flow of the effective coupling constants is 
not diverging in the infrared region. We stress that our 
method applies to a large class of perturbations of the 2d 
Ising model, and for sake of definiteness we restrict our 
analysis to AT. 

In order to present our result, we find convenient to 
introduce the variables 



t = 



id) +t(2) 



(2) 



with t'^^^ — tanhj'^^^ j = 1,2. The parameter t has 
the role of a reduced temperature and u measures the 
anisotropy of the system. We shall consider the free en- 
ergy or the specific heat as functions of t,u,l. When 
I = the specific heat Cy can be immediately computed 
from the Ising model exact solution; Cv is diverging at 
t = tf = \f2 — 1 ± |it| and near the the critical temper- 
atures the specific heat shows a logarithmic divergence: 
Ct, ~ — Clog|t — i^l, where C > 0. If the anisotropy 
is strong the two Ising subsystems have very different 
critical temperatures, hence one can expect that if one 
system is almost critical the second one will be out of 
criticality; then mean field arguments based on the fact 
that two Ising are coupled by a density-density interac- 
tion suggest that the effect of the coupling is to change 
at most the value of the critical temperatures. On the 
other hand if the anisotropy is small the two system will 
become critical almost at the same temperature and the 
properties of the system could change drastically. 

Our main result is the following theorem; the detailed 
proof can be found in 0, |23 . 



Theorem. For I small enough the AT model admits 
two critical points of the form: 

tfil, u)^V2-l + v{l) ± l^l^+^ll + 5{l, u)) . 



Here v and b are 0(V) corrections and rj = —bl + 0(P) 
with b > 0. If t ^ tf the free energy of the model is 
analytic in l,t,u and the specific heat Cv is equal to: 



-FiA^"^" log 



\t-t-\-\t-ti 
A2 



1 - A^'fc 



+ F3, (3) 



where: 2IS? = (t - t^)^ + (t - t+ f; r]c = al + 0{P), 
a 7^ 0; and Fi, F2, F3 are functions of t,u,l, bounded 
above and below by 0(1) constants. 

1) First note that the location of the critical points is dra- 
matically changed by the interaction. The difference of 
the interacting critical temperatures normalized with the 
free one G(/, u) = {t+ {I, u)-t- {I, u))/{t+ (0, u)-t- (0, u)) 
rescales with the anisotropy parameter as a power law 
~ \u\^ , and in the limit u — s- it vanishes or diverges, de- 
pending on the sign of I (this is because 77 = —bl + 0{P), 
with 6 > 0). In Fig. 1 we plot the qualitative behaviour 
of G{1, u) as a function of it, for two different values of I 
(i.e. we plot the function u^, with 77 = 0.3,-0.3 respec- 
tively) . 




FIG. 1: The behaviour of the difference G between the in- 
teracting critical temperatures normalized to the free one, for 
two different values of I; depending on the sign of the inter- 
action, it diverges or vanishes in the isotropic limit. 

As far as we know, the existence of the critical index 
rjil) was not known in the literature, even at a heuristic 
level. 

2) There is universality for the specific heat, in the sense 
that it diverges logarithmically at the critical points, as 
in the Ising model. However the coefficient of the log 
is anomalous: in fact if t is near to one of the critical 
temperatures A ~ \/2|m|^+'' so that the coefficient in 
front of the logarithm behaves like ~ with 
r]c a new anomalous exponent 0{l); in particular it is 
vanishing or diverging as u ^ depending on the sign of 
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I. We can say that the system shows an anomalous uni- 
versality which is a sort a new paradigmatic behaviour: 
the singularity at the critical points is described in terms 
of universal critical indexes nevetheless in the isotropic 
limit u — *■ 0, some quantities, like the difference of the 
critical temperatures and the constant in front of the 
logarithm in the specific heat, scale with anomalous 
critical indexes, and they vanish or diverge, depending 
on the sign of I. 

3) Eq(3) clarifies how the universality-nonuniversality 
crossover is realized as u ^ 0. When u ^ only the 
first term in eq(3) can be log-singular in correspondence 
of the two critical points; however the logarithmic term 
dominates on the second one only if t varies inside an 
extremely small region 0(|u|^+''e~'^/l'l) around the crit- 
ical points (here c is a positive 0(1) constant). Outside 
such region the power law behaviour corresponding to 
the second addend dominates. When u — > one recovers 
the power law decay found in the isotropic case 



Cv — -Fs 



1 - |t - teP"- 
Vc 



In Fig. 2 we plot the qualitative behaviour of Cv as a 
function of t. The three graphs are plots of eq(3), with 
Fi = Fa = 1, F3 = 0, M = 0.01, ry = r/c = 0.1,0,-0.1 
respectively; the central curve corresponds to the case 
= 0, the upper one to rj < and the lower to 77 > 0. 
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FIG. 2; The behaviour of the specific heat Cv for three dif- 
ferent values of I, showing the log-singularities at the critical 
points; in the isotropic limit the two critical points tend to- 
coincide, the lower curve becomes continuous while the upper 
develops a power law divergence. 

We now sketch the proof of the above Theorem (for 
a detailed proof we refer to [H IH). We start from 
the well known representation of the Ising model free en- 
ergy in terms of a sum of Pfa ffian s | 2 4l which can be 
equivalently written (see Ref. |l7l as Grassmann 



functional integrals, formally describing massive non in- 
teracting Majorana fcrmions ?/;, 1/; on a lattice with action 

\ [Mdi - »5o)Vx + ^x(9i + - (4) 

X 

-2ili^^{di + ao)^x] + i{V2 ~ 1 - <)V^^Vx , 

where dj are discrete derivatives; criticality corresponds 
to the massless case. If Z = the free energy and specific 
heat of the AT model can be written as sum of Grass- 
mann integrals describing two kinds of Majorana fields, 
with masses m^^^ = t^^)- VS-hl and m^^^ = t^'^'^ - \/2+l. 

If Z ^ again the free energy and the specific heat can 
be written as Grassmann integrals, but the Majorana 
fields are interacting with a short range potential. By 
performing a suitable change of variables and integrating 
out the ultraviolet degrees of freedom, the effective action 
can can be written as 

x,0,a 

+z0A*iV'0,xV'-^,_x + ^iV'txV'i:xV'ii,xV'"i,x] + Wi 

where a = ± is a creation-annihilation index and = ±1 
is a quasi-particle index, ui and have the role of 
two masses and it holds cri = 0{t ~ + I) + 0{l), 
fii = 0{u). Wi is a sum of monomials of ■0 of arbitrary 
order, with kernels which are analytic functions of Zi; an- 
alyticity is a very nontrivial property obtained exploit- 
ing anticommutativity properties of Grassman variables 
via Gram inequality for determinants. The ip'^ are Dirac 
fields, which are combinations of the Majorana variables 

j = 1, 2, associated with the two Ising subsys- 
tems. 

One can compute the partition function by expanding 
the exponential of the action in Taylor series in I and 
naively integrating term by term the Grassmann mono- 
mials, using the Wick rule; however such a procedure 
gives poor bounds for the coefficients of this series that, 
in the thermodynamic limit, can converge only far from 
the critical points. 

In order to study the critical behaviour of the system 
we perform a multiscale analysis involving non trivial re- 
summations of the perturbative series. The first step is to 
decompose the propagator ^(k) as a sum of propagators 
more and more singular in the infrared region, labeled 
by an integer h < 1, so that g{k) ^ EL-oo 3*''^^), 
^('')(k) ~ 7^'*. We compute the Grassmann integrals 
defining the partition function by iteratively integrat- 
ing the propagators g^^\g''^\ ■ ■ ■ After each integration 
step we rewrite the partition function in a way simi- 
lar to the last equation, with Z^, cr/j, /i/j, Z/i, Wh replacing 
Zi, CTi, /ii, Zi, Wi, in particular the masses and the wave 
function renormalization are modified; the structure of 
the action is preserved because of symmetry properties; 
moreover Wh is shown to be a sum of monomials of 
of arbitrary order, with kernels decaying in real space on 
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scale 7"^^, which are analytic functions of {l^, . . . , h}, if 
Zfe are small enough, k > h, and |crfc|7~'^, |/ifc|7~''' < 1; 
again analyticity follows from Gram-Hadamard type of 
bounds. 

All the above construction is based on the crucial prop- 
erty that the effective interaction at each scale does not 
increase \lh\ < 2|/|; such property is a consequence of 
the validity of some non perturbative approximate Ward 
identities [23; "approximate" refers to the fact that, be- 
cause of the presence of masses and of an ultraviolet cut- 
off, the Ward identities are different from the usual for- 
mal ones; the error terms are shown to be small, in a 
suitable sense. For ah, fJ-h, Zh, we find that, under the 
iterations, they evolve as: at — aij''^"'-, — ^17"''^"', 
Zh — 7~^^' with 61, 62 explicitely computable in terms 
of a convergent power series. 

We perform the iterative integration described above 
up to a scale /i^ such that (|(T/i* | + |)7~''i = 0(1). For 
scales lower than we return to the description in terms 
of the original Majorana fermions ^/;(2,<'!i) as- 

sociated with the two Ising subsystems. One of the two 
fields (say is massive on scale hi (so that the 

Ising subsystem with j = 1 is "far from criticality" on the 
same scale); then we can integrate the massive Majorana 
field without any further multiscale analysis, ob- 

taining an effective theory of a single Majorana field with 
mass |cr/ij| — 1/^/1*1, which can be arbitrarly small; this is 
equivalent to say that on scale hi we have an effective de- 
scription of the system as a single perturbed Ising model 
with anomalous parameters near criticality. The inte- 
gration of the scales < hi is performed again by a mul- 
tiscale decomposition similar to the one just described; 
an important feature is however that there are no more 
quartic marginal terms, because the anticommutativity 
of Grassmann variables forbids local quartic monomials 



of a single Majorana fermion. Criticality is found when 
the effective mass on scale —00 is vanishing; the values 
oi t,u for which this happens are found by solving a non 
trivial implicit function problem. 

Technically it is an interesting feature of this problem 
that there are two regimes in which the system must be 
described in terms of different fields: a first one in which 
the natural variables are Dirac Grassmann variables, and 
a second one in which they are Majorana; the scale hi 
separating the two regimes is dynamically generated by 
the iterations. In the first regime the two entangled Ising 
subsystems are undistinguishable, the natural description 
is in terms of Dirac variables and the effective interaction 
is marginal; in the integration of such scales nonuniver- 
sal indexes appear. In the second region the two Ising 
subsystems really look different, one appear to be (al- 
most) at criticality and the other far from criticality on 
the same scale; the parameters of the two subsystems 
are deeply changed (in an anomalous way) by the previ- 
ous integration; in this region the effective interaction is 
irrelevant. 

In conclusion we have presented some new rigorous 
results on the critical behaviour in the Ashkin-Teller 
model, for weak coupling and any value of the anisotropy. 
Via multiscale integration methods we have computed 
the specific heat and the location of the critical temper- 
atures in terms of convergent power series and we have 
predicted the existence of an unknown critical exponent 
describing the scaling of the gap between the critical 
temperatures in the isotropic limit. Moreover we gave 
a detailed description of the crossover between the uni- 
versal critical behaviour holding in the anisotropic case 
and the anomalous nonunivcrsal behaviour holding in the 
isotropic limit. 
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